Non collinear magnetism and single ion anisotropy in multiferroic perovskites 
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The link between the crystal distortions of the perovskite structure and the magnetic exchange 
interaction, the single-ion anisotropy (SIA) and the Dzyaloshinsky-Moriya (DM) interaction are 
investigated by means of density-functional calculations. Using BiFeOs and LaFeOs as model sys- 
tems, we quantify the relationship between the oxygen octahedra rotations, the ferroelectricity and 
the weak ferromagnetism (wFM). We recover the fact that the wFM is due to the DM interaction 
induced by the oxygen octahedra rotations. We find a simple relationship between the wFM, the 
oxygen rotation amplitude and the ratio between the DM vector and the exchange parameter such 
as the wFM increases with the oxygen octahedra rotation when the SIA does not compete with 
the DM forces induced on the spins. Unexpectedly, we also find that, in spite of the d 5 electronic 
configuration of Fe 3+ , the SIA is very large in some structures and is surprisingly strongly sensitive 
to the chemistry of the A-site cation of the ABO3 perovskite. In the ground R3c state phase we 
show that the SIA shape induced by the ferroelectricity and the oxygen octahedra rotations are in 
competition such as it is possible to tune the wFM "on" and "off" through the relative size of the 
two types of distortion. 



I. INTRODUCTION 

During the last ten years, there has been a huge in- 
crease of interest in developing magnetoelectric multifer- 
roic materials. In such materials, ferroelectric and mag- 
netic ordering coexist together and can be coupled such 
that the magnetization is affected by an electric field and 
the polarization by a magnetic field. The reason for this 
interest is related to impact potential in technological 
applications in transducers, attenuators, filters, informa- 
tion storage and spintronics 1 5 . The ideal magnetoelec- 
tric multiferroic would be a compound in which a large 
spontaneous polarization were coupled with a large mag- 
netization so that flipping the former could flip the lat- 
ter and vice- versa. However, such an ideal compound is 
not known today. A particularly promising direction is 
to switch magnetization by 180° using an electric field 
in materials exhibiting weak ferromagnetism (wFM) 6 ' 7 . 
In weak ferromagnets, the magnetization is small, but it 
has been proposed theoretically that in the presence of 
a ferroelectric polarization one can switch the wFM by 
reversing the polarization 6-9 . 

Following this pathway for magnetoelectric control of 
wFM by electric polarization, it is crucial to understand 
the underlying mechanisms that link the wFM to the 
crystal distortions. Dzyaloshinsky and Moriya (DM) 
showed that spin-orbit interaction (SOI) mediates a spin- 
spin coupling of the form D.S\ x 62 (the so-called DM 
interaction) that is usually responsible for wFM 10 ' 11 . 
Moreover, Bertaut showed by symmetry considerations 
that under sufficiently low crystal symmetry, SOI also 
permits the single ion anisotropy (SIA) to cause non- 
collinear magnetic arrangements 12 . However, the latter 
interaction has been less well studied, and in most cases 
only the DM interaction has been considered to be re- 
sponsible for the wFM. 

In multiferroic perovskites, the initial para-electric cu- 
bic reference structure can be deformed by the presence 



of two main lattice instabilities: the antiferrodistortive 
(AFD) instabilities, consisting of non-polar oxygen oc- 
tahedra rotations, and the ferroelectric (FE) instability, 
responsible for the polarisation. It has been recognized 
that spin canting is induced by the AFD distortions and 
thus the wFM is directly linked to the amplitude of the 
AFDs 13,14 . From symmetry analysis, some of such dis- 
torted structures do not allow for spin canting while some 
allow only the DM to be responsible for the spin cant- 
ing and a few cases allow spin canting through both DM 
and SIA 12 . In spite of this, in multiferroic perovskites, 
no systematic study has been performed to analyze the 
details of the coupling between spin canting and lattice 
distortions. 

The aim of the present study is to quantify from first- 
principles calculations the mechanisms leading to wFM 
in two representative perovskites, LaFeOs (LFO) and 
BiFe03 (BFO). Both compounds show AFD distortions 
(of different types), and BFO also has a FE instability, 
whereas LFO does not. To understand the links between 
these lattice distortions and the magnetism, we decom- 
pose the magnetic interaction into three main types, ex- 
change, DM and SIA and look at how they are affected 
by the amplitude and combinations of the different lat- 
tice distortions. This systematic analysis allows us to 
understand the coupling between structural distortions 
and spin canting, and propose some guidelines for the 
design of magnetoelectric wFMs through the SIA. 



II. TECHNICAL DETAILS 

The first-principles calculations were performed us- 
ing density-functional theory (DFT) as implemented in 
the VASP code within the Projector Augmented Wave 
(PAW) method 15 ' 16 . Local spin density approximation 
with an additional Hubbard (LSDA+C/) was used for 
the exchange-correlation functional. The Hubbard pa- 
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rameter U and the exchange interaction J that treat the 
Fe d electrons were set to U=5 eV and J=0 eV, values 
that have been shown to be optimal for Fe 3+ in LaFeOa 
and BiFeOs 13 ' 17,18 . The wave functions were expanded 
in plane waves up to a cut-off energy of 500 eV. For in- 
tegration of the Brillouin-zone of the supercell made of 
2x2x2 cubic perovskite units, we used a 3 x 3 x 3 
Monkhorst-Pack k-point mesh shifted by \ x | x | . The 
convergences were tested with a cut-off energy up to 700 
eV and a 6 x 6 x 6 k-point mesh, with no significant 
change on the calculated SIA. Spin-orbit coupling (SOC) 
was included in calculating both total energies and forces 
to incorporate the coupling between the spins and the 
lattice. Symmetrization was switched off to remove any 
artificial constraints on the possible spin ordering. 

AFD distortions are described by rotating the oxygen 
octahedra around the central Fe 3+ ion, where the rota- 
tion axis is one of the Cartesian axis. We use the modi- 
fied Glazer notation 19 , where the complete AFD pattern 
is classified by a triplet a a 6^c 7 with a, b and c refering 
to the amplitude of rotation around the Cartesian x, y 
and z-axis and the superscripts a, /3,7 = {+, —,0} refer 
to the type of rotation, + and — stand for in-phase and 
out-of-phase respectively and for no rotation in this 
direction. 

To control the amplitude of AFD rotations, we treated 
the Fe-0 bond length as constant when freezing in the 
oxygen octahedra rotations, which is a good approxima- 
tion with respect to the fully relaxed AFD structures. 
The rotation is sometimes thought of as the displacement 
of the oxygen along the edge of the cubic unit cell, which 
is pictured by the red arrows in Fig.l.b. However, this 
view is valid for very small oxygen octahedra rotation 
only since for large amplitude of rotations, if the oxygen 
are kept on the edge of the cubic cell, this will strongly 
stretch the Fe-0 bonds in the plane of rotation, a dis- 
tortion that is not observed in fully relaxed structures. 
Instead, if the Fe-0 bond lengths are kept constant (we 
chose 3.90 A, which is close to the Fe-0 bond lenght 
observed in the ground state of BFO 20 and LFO 21 ), an 
oxygen octahedra rotation is accompanied by a geometric 
shrinkage of the cell parameter in the plane perpendic- 
ular to the rotation axis as illustrated in Fig.l.b and c. 
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FIG. 1. Schematic view of the geometrical skrink of the per- 
ovskite cell in the presence of an AFD distortion along z-axis 
perpendicular to the paper plane. 

The interaction parameters, exchange J^, DM D^, 



and SIA <1>^, in the Hamiltonian 12 

H = -2[JijSi • Sj + Dij • (Si x Sj) +Si-$i- Si] 

were extracted from first-principles calculations. We used 
a 2 x 2 x 2 supercell with 8 ABO 3 formula units. To sep- 
arate these three interactions, we used two approaches. 
First, we performed artificial calculations where selected 
Fe 3+ ions were replaced by nonmagnetic Al 3+ ions so 
that only one of the three interactions was retained. For 
example, the SIA is extracted by replacing all except one 
of the Fe cations with neutral Al cations and then per- 
forming constrained calculations of the direction of the 
remaining spin (Lagrange multiplier) to resolve the en- 
ergy surface. The data points were then fitted to extract 
the parameters Ki on the usual expressions for the SIA 22 . 
By removing all the surrounding magnetic cations, we as- 
sure the vanishing of all magnetic interactions (exchange 
and DM) except for the SIA 23 . We note that in this spe- 
cial configuration the point symmetry of the remaining 
Fe atom is the same as in the configuration where all the 
Fe atoms are present, a condition that guarantees that 
the crystal field splitting will be similar in the two config- 
urations. The exchange parameters (Jj) and DM vectors 
(Dij) are extracted by replacing all except for two of 
the Fe cations with Al. The exchange constant can then 
be extracted directly from the energy difference between 
parallel and anti-parallel arrangements of the two spins 
using 

a=x,y,z 

where S is the spin moment (//#), E[s a ,s a ] is the total 
energy of a spin configuration and s a the spin direction. 
The DM vectors can be extracted by perpendicular ar- 
rangements of the two spins using 

[Di j } a = ^(E[§isi]-E[§i,§i}) 

with a, b and c being three perpendicular orientations. 
This method is valid only if both magnetic ions have the 
same non-uniaxial SIA or if the SIA is uniaxial, which is 
the case in the structures analysed in the present study. 
A second approach is that proposed by Xiang et al, 24 
which does not have this restriction on the symmetry 
of the SIA, but requires twice as many calculations. We 
compared the results from both methods and find simular 
results. 



III. EXCHANGE INTERACTION 

In Tab. I we report the calculated exchange parameters 
Jij extracted from our DFT calculations for LaFeOs and 
BiFeOs in cubic, Prima, R3c and R3c structures (cell 
parameters and amplitude of distortions are kept fixed 
to the same values for LFO and BFO in each case). The 



3 



cubic structure always shows stronger exchange interac- 
tions than the distorted structures because of the large 
orbital overlap associated with its 180° Fe-O-Fe bonds. 
It is interesting to see that the J's for BFO and LFO are 
very close in all the structures, we note, however, that in 
BFO the exchange parameters are always smaller than 
in LFO due to the fact that both BFO and LFO have 
the same Fe 3+ magnetic cation and we are imposing the 
same environment. We also remark that the J's decrease 
with increasing the amount of either AFD or FE struc- 
tural distortion. In addition to the method described in 
the previous chapter, we also computed these exchange 
parameters with the method proposed by Xiang 24 which 
gives similar results. For example we find with this sec- 
ond method J ac = 6.70 meV and 6.06 meV respectively 
for the Prima and R3c phase of LFO. 
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cubic 


0°0°0° 
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7"8 + 7" 


6.74 [6.70] 


6.83 




R3c 


9"9"9" 


6.17 [6.06] 






R3c 


9-9-9-+1.0 FE 


5.64 




BFO 


cubic 


0°0°0° 


7.36 






Prima 


7-8+7- 


6.52 


6.68 




R3c 


9-9-9- 


5.96 






R3c 


9-9-9-+I.O FE 


5.36 





TABLE I. Our calculated Heisenberg exchange constants 
(meV) for LaFeOs and BiFeOs in different crystal structures. 
We used a modified Glazer notation to indicate the amplitude 
of the oxygen octahedra rotations (in degree) in the three di- 
rections. The amplitude of the FE distortions in the R3c 
structures (1.0 FE) are those given by the minimum of en- 
ergy when freezing in the polar unstable mode of the corre- 
sponding paraelectric reference R3c structure. The numbers 
in brackets were obtained with the method of Xiang et al. 24 . 

Because all the exchange interaction constants are pos- 
itive, the favored spin structure is antiferromagnetic G- 
type for all the structures. This is confirmed by looking 
at the energy differences between the G, G, A and F 
types of magnetic orders reported in Tab. II, where it is 
clear that the lowest energy state is always G-type AFM 
order. 



Magnetic 
ordering 


cubic 


LFO 

Puma 


cubic 


BFO 

Puma 


A-type 


-0.12 


-0.10 


-0.11 


-0.10 


C-type 


-0.22 


-0.19 


-0.21 


-0.19 


G-type 


-0.31 


-0.27 


-0.30 


-0.27 



TABLE II. Energy difference per formula unit (in eV) of the 
three main antiferromagnetic orders with respect to the fer- 
romagnetic F-type configuration. 



Next we calculate how the DM interaction is affected 
by the same structural distortions. 



IV. DM INTERACTION 

In this section we compare the amplitude of the DM 
vectors between BFO and LFO in Prima, R3c and R3c 
structures. 



A. Prima structure 

The symmetry elements of the orthorhombic Prima 
space group and the antisymmetry of the DM interac- 
tion (Dij = —Dji) determine that the DM vectors have 
the form shown in Fig. 2 described by five parameters: 
<^ac, /^ac, 7oc, <^6, 76 - 24 - Note that the Dy vectors corre- 
sponding to two Fe moments are located on the inter- 
mediate oxygen ions. The DM vector at each O site 
is known to be perpendicular to the mirror plane run- 
ning through the associated Fe-O-Fe unit 10 , and given 
by oc Xi x x j, where Xk is the bonding vector 0-Fe& 




FIG. 2. DM vectors located at the oxygen atom between two 
Fe for the Puma structure. 

In Tab. Ill we report our calculated parameters of the 
DM vectors for both LFO and BFO. The absolute values 
of the components of the total DM vector, which are the 
sums of the local DM vectors shown in Fig. 2, are reported 
in Tab. III. The amplitude of the DM vectors are of the 
order of several hundred /ieV. Comparing the values with 
those calculated for LFO by Kim et al. 14 , we find that 
our calculated absolute values are smaller but the ratio 
Dac/Jac = 0.018 and D^/J^ = 0.017 are consistent with 
the values of Kim et al. who reported 0.020 and 0.021 
respectively. One possible reason for the difference could 
be the different value of U used (7.5 eV in the microscopic 
model of Ref.14). In Fig. 3 we show the corresponding di- 
rections for these calculated DM vectors in LFO (blue 
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FIG. 3. DM vectors of two adjacent Fe ions along the b 
direction shown from two different angles. The directions of 
the geometrically constructed (calculated) DM vectors Dij 
are shown in red (blue). 



arrows) which are in good agreement with those we ex- 
pect from geometrical constructions (red arrows). 
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TABLE III. DM vector components (/xeV) of LFO and BFO 
in the Prima phase (7~8 + 7~), calculated in this work (top 
two rows) and in 14 (bottom row). 



The DM interaction between two spins Si and sj can 
be seen as a force induced by the spin Sj on the spin 
Si and can be expressed as fV = Dy x Sj. 25 The to- 
tal force acting on the spin Si due to its six neighbor- 
ing spins sj is thus fi = Y^=i^\- Using this force 
analysis of the DM interaction, in the Prima struc- 
ture with G 2 -type ordering (G-type where the spins lie 
along the z direction), one finds two types of force hav- 
ing the same norm: fi = (— 4/3 ac , — 2a^ — 4a ac , 0) and 
f2 = (±4/3 ac , — 2a& — 4a ac ,0). These two types of force 
differ only in their orientation along the x direction and 
show that half of the spins feel a force along +x and 
half along —x. Therefore we see that the DM interac- 
tion parameter f3 ac causes the A-type canted AFM pre- 
dicted by symmetry along the x direction (A x ). Along 
the y direction all the spin sites feel a force in the same 
direction which causes the wFM canting along the y di- 
rection (F y ) with a strength determined by a& and a ac . 
Along the easy z direction, no force is induced by the 
DM interaction. This is in good agreement with simple 
symmetry considerations that show that G z , A x and F y 
have the same symmetry transformation in the Prima 
perovskite structure 12,26,27 . The resulting torques are 
(±364, -446, 0) fieV for LFO and (±208, -454, 0) /xeV 
for BFO. The forces induced by the DM interaction are 
of the same amplitude along the y direction [F y canting) 
in both LFO and BFO, while the resulting force along 
the x direction (A x canting) is larger in the case of LFO 



than in BFO. 



B. DM in R3c and R3c 

In R3c and i?3c, the symmetry of the crystal struc- 
ture requires that the total DM vector is along the [111] 
direction 13 . The direction of canting is related to the 
sign of D in such a way that the three vectors D, Si 
and S2 build up a right-handed system. The three lo- 
cal DM vectors (Fig. 4) can be fully described using only 
two independent parameters a and /?: Di = ((3, a, a), 
D 2 = (a,/3,a), D 3 = (a,a,/3). 

In Tab. IV we report our calculated values of the a 
and j3 parameters and the total DM vector for LFO and 
BFO in both R3c and R3c phases, a is always much 
larger than (3 in LFO while in BFO they have the same 
amplitude. The total magnitude of the DM vectors are 
however similar for LFO and BFO as well as in the two 
structures R3c and R3c. The FE distortions have the 
tendency to reduce the DM vector but with a relatively 
small effect. 
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BFO 


52 


66 


170 


R3c 


LFO 


50 


20 


120 




BFO 


48 


50 


146 



TABLE IV. Calculated a and /3 parameters (/xeV) and mag- 
nitude of the DM vector components (/xeV) of LFO and BFO 
in the R3c and R3c phases. 



As we did for the Prima phase, we can compute the 
forces induced on the spins due to the DM interaction. 
In the case of R3c one finds that all the spin sublattices 
feel the same force f = (— j3 — 2a, —/3 — 2a, j3 ± 2a). This 
force is perpendicular to the spins and lies in the (111) 
plane which causes a canting of all the spins in the (111) 
plane. These results are in good agreement with the ex- 
perimental data where the wFM moment is found to be 
in the plane perpendicular to the [111] direction. 28,29 




FIG. 4. Sketch of the orientation of the three local DM vectors 
in R3c and R3c structures. 
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In the Prima and R3c structures, we find that the DM 
interaction has the tendency to cant the spins away from 
the easy axis, related to the wFM or wAFM observed 
in these structures. In the next section we calculate the 
effect of the SIA on the magnetic ground state of BFO 
and LFO. 



V. SINGLE ION ANISOTROPY 

Finally, we calculate the behaviour of the SIA energy 
in the cubic and different distorted structures, focusing 
in particular on how the AFD and FE distortions and 
their combinations affect SIA. Because the SIA has a 
more complex link with the crystal distortions than the 
DM and exchange interactions, we analyse the SIA for 
a°a°c + , a°a°c _ , a°6 + 6 + , a°b~b~, in addition to cubic, 
Prima, R3c and R3c structures. 

A. Cubic perovskite 

A perovskite without any oxygen octahedra rotation 
has the cubic space group Pm3m 19 . For cubic crystal 
field splitting, the SIA can be described by 4 th and 6 th 
order terms: 

E^tiou) = Kt{ala 2 y + a 2 y a 2 z + a 2 z a 2 x ) + K c 2 (a 2 x a 2 y a 2 z ) 

(1) 

where is the normalized projection of the spin in the 
i-direction with the constraint Y^=i a l = 1 anc ^ an d 
K 2 are the SIA parameters. Our calculated SIA con- 
stants obtained for LFO and BFO in the cubic structure 
are given in Tab.V. 





Kl 


K% 


LFO 


-1.62 


0.00 


BFO 


-3.66 


0.06 



TABLE V. Calculated SIA constants of LFO and BFO in the 
cubic structure fitted to Eq.l. 



As we can see, K° is negative for both BFO and LFO, 
which indicates that the spins point in any of the three 
diagonal directions [111] 22 . As expected for cubic sym- 
metry, the anisotropy energy is very small (a few jaeV). 
Interestingly, it is two times larger in BFO than in LFO. 
The 6 th order anisotropy constant K% is close to zero for 
LFO and is two orders of magnitude smaller than K\ in 
BFO. This shows that, in the cubic structure, the chem- 
istry of the A-site ion affects the anisotropy at the B site, 
with Bi causing a larger anisotropy than La. 

B. a°a°c + and a°a°c~ structures 

In this section, we consider the effect of a single in- 
phase AFD (a°a°c + ) and out-of-phase AFD (a°a°c~) on 



the SIA energy. Since the symmetry is tetragonal, we use 
the following general expression to fit our SIA energy 22 : 

Esia {P,<j>)= K x sin 2 (0) + K[ sin 2 (0) cos(20) (2) 
+K 2 sin 4 (0) + if 2 sin 4 (0) cos(40) (3) 

where G [0, tt] is the polar angle between the spin direc- 
tion and the local z-axis and (j) G [0, 2tt] is the azimuthal 
angle in the plane perpendicular to = 0. 

In Tab. VI we report the values of the SIA constants 
for AFD rotations of 10° about the z axis (0°0°10 + and 
0°0°10") in both LFO and BFO. In the two compounds 
we obtain K\ » K 2 > K 2 > for a°a°c+, while 
for a°a°c~ K\ > only for BFO. According to Eq.3, 
Ki > means that the anisotropy is predominantly uni- 
axial. Because of the small size of K 2 and K'i relative 
to Ki, a second-order uniaxial model serves as a good 
description for the SIA. It is interesting to see that the 
larger anisotropy of BFO versus LFO reported in the cu- 
bic structure is further emphasized in the presence of an 
oxygen octahedral rotation: K\ is several hundred fieV 
for BFO while it is only few /ieV for LFO. The uniax- 
ial anisotropy of BFO is very robust against the oxygen 
AFD distortions, however we remark that in LFO a°a°c + 
gives rise to uniaxial SIA [K\ > 0) while a°a°c~ gives 
rise to easy plane SIA {K\ < 0). This can be related 
to the small value of the SIA of LFO that can be easily 
affected by small structural changes. 
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LFO 


6.7 
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BFO 


264.0 





3.5 
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0°0°10" 


LFO 


-1.6 





1.3 


0.2 




BFO 


235.3 





4.3 


0.8 



TABLE VI. SIA constants (/ieV) for LFO and BFO with 
a°a° b + and a°a°b- types of rotation (rotation amplitude of 
10°). 



In Fig. 5. a we plot the energy variation when the spin 
is turned in the yz plane (xz plane being equivalent) of 
LFO. As we can see, the lowest energy is obtained for 
spins orientated along the z axis, which means that when 
one AFD rotation develops in the perovskite structure, 
the rotational axis becomes the easy axis for the spins. 
The SIA in the 7/z-plane (or xz-plane) with the presence 
of an AFD rotation along the z axis is larger than that 
show for the perfect cubic cell but still stays small in LFO 
(6.7 fieV) while it is sizeable in BFO (264 /ieV). 

In Fig.5.b we also show the SIA energy of LFO in 
the plane perpendicular to the oxygen rotation axis (xy 
plane). One can see the spin does not point in the same 
direction on Fe atom site A (where the oxygen rotation 
is counter clockwise around Fe) and B (where the oxy- 
gen rotation is clockwise around Fe). On atom site A 
the energy minimum is at 35° (blue squares on Fig.5.b) 
while for atom site B the energy minimum is at 55° (red 
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FIG. 5. SI A energy of LFO in (a) the zy-plane of LFO for 
a°a°c + (blue squares) and a°a°c~ (red circles), (b) in the xy 
plane in the presence of a°a°c + (blue squares for SI A at site 
A and red circles for site B, see text) and (c) SI A energy of 
BFO in the zy plane in the presence of a°a°c + (blue squares) 
and a°a°c- (red circles) . In panel (b) we also report the total 
SIA energy (green stars) made of the sum of all the local SIA 
of each Fe sites. We also compare in panel (b) the magne- 
tocrystalline energy calculated by turning all the spins (no Al 
replacement) in the xy plane with the same angle. In panels 
(a) and (c), the angle 0° represents the [001] direction (z) and 
the angle 90° represents the [010] (y) direction. In panel (b) 
the angle 0° represents the [100] direction and the 90° an- 
gle represents the [010] direction. Note the energy amplitude 
differences between panels (a), (b) and (c). 



diamonds on Fig.5.b). This clearly show that the local 
spin anisotropy and thus the local spin direction follows 
directly the AFD rotation amplitude (45° ±10°). How- 
ever, the global anisotropy is given by the sum of the 
SIA of all of the magnetic cation sites. This sum of SIAs 
is represented in Fig.5.b by the purple triangles. The 
global anisotropy gives rise to a minimum of energy at 
45°, highlighting the fact that the shift of ±10° with re- 



spect to the [110] direction of each site compensate each 
other such that the global anisotropy is lowest in the [110] 
direction. In Fig.5.b we also report the global SIA calcu- 
lated by turning simultaneously the spins with all A and 
B site occupied with Fe atoms (green dots). With all the 
spins together, we recover the energy minimum at 45° 
and with energy amplitude in good agreement with the 
sum of the single spins result, highlighting the correct- 
ness of our approximation to replace all the surrounding 
Fe atoms by Al atoms. It is also interesting to see that, 
apart from a phase shift of 10°, in the plane xy perpen- 
dicular to the AFD rotation axis, the amplitude of the 
anisotropy is similar to that found in the cubic structure. 
We note that this global SIA is the magnetocrystalline 
anisotropy (MCA). Comparing with BFO, the SIA en- 
ergy in the zy plane (Fig.5.c, red triangles) gives the same 
easy axis (the [001] direction minimizes the energy) but 
with a much larger amplitude (close to 300 fieV). 

Performing the same SIA analysis with the a°aPc~ 
structure and the same angle of 10° of oxygen octahe- 
dral rotation, we recover similar results as in the a°a°c + 
case but with some differences between BFO and LFO. 
In BFO, the shape and amplitude of the SIA in a°a°c~ 
is the same as in a°a°c + , a result not surprising since 
locally, for one single spin, the first nearest-neighbors are 
the same in both a°a°c~ and a°a°c + structures: In the 
xy-pl&ne we have the same alternating A and B sites 
while in the z direction the A and B sites, are inverted 
in a°a°c~ (out-of-phase rotation along the z direction) 
with respect to a°a°c + (in-phase rotation along the z di- 
rection). However this does not hold for LFO since the 
easy axis changes from [001] in a°a°c + to [111] in a°a°c~ . 
This can be related to the fact that the anisotropy is 
much smaller in the case of LFO (few /xeV) and so small 
changes in the structure, which give rise to changes in 
the SIA of a few /ieV, can change the easy axis. 



C. oPb + b + and a°b b structures 

In this section we analyze the effect of two in-phase and 
two out-of-phase AFD rotations (a°6 + 6 + and a°b~b~) on 
the SIA. We consider angles of AFD rotation of the same 
amplitude (0°10+10+ and 10 _ 10-) and we compare 
with the previous case where we considered only one ro- 
tation. The structures of a°6 + 6 + and a°b~b~ have space 
groups 14/ m and Imma respectively and the model that 
describes the SIA is that given by Eq.3. 

In Tab. VII we report the fitted K parameters obtained 
in the 0°10 + 10+ and 0°10"10" structure. The quanti- 
zation axis (0 = <p = in Eq.3) is the local hard axis, 
that is the [011] direction for all the cases reported in 
Tab. VII. In both compounds we obtain an easy plane SIA 
(Ki is negative). As in the a a %+/- 

cases, the second- 
order anisotropy constants K\ and K[ are larger than the 
fourth-order K2 and K^. Again, we find the anisotropy 
of BFO (474.7 /ieV) to be much higher than that of LFO 
(130.2 /ieV) for both, out-of-phase and in-phase AFD dis- 
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tortions. Interestingly, and in contrast to the case with 
oxygen rotation in only one direction, the SIA is strongly 
non-uniaxial since it is modulated by the large value of 
the K[ constant. We also remark that the SIA energies 
increase when adding a second AFD distortion to the 
a°a°6 + / - structures. 









K x 


K[ 


K 2 


K'i 


0°10+ 


10+ 


LFO 


-71.8 


24.1 


-1.6 









BFO 


-515.6 


201.5 


-2.9 





0°1(T 


10" 


LFO 


-130.2 


27.2 


-1.4 









BFO 


-474.7 


204.9 


-2.7 






TABLE VII. Fitted SIA constants (ueV) of LFO and BFO in 
a 6+ 6+ and a°b~b~ structures with 6=10°. 



D. a b + a Pnma structure 

As for the previously reported distorted perovskites, 
the SIA of the Pnma phase is also described by Eq.3. Our 
calculated coefficients for LFO and BFO in their Pnma 
phase are given in Tab. VIII. As obtained for the previous 
cases with one (a°a°c+ and a°a°c~) and two (a°6+6+ and 
a°b~b~) oxygen rotations, the K\ constant is the largest. 
However, and in contrast to these previous strucures, in 
the Pnma phase K\ is negative, which means that the 
shape of the SIA is mainly of an easy plane rather than an 
easy axis type. Here again we remark that the anisotropy 
of BFO is much larger than that of LFO. For LFO K[ is 
small (7/ieV) while it is larger for BFO (32 /ieV), showing 
that the SIA for both LFO and BFO slightly deviate from 
a perfect easy plane. The K 1 ^ and are zero for BFO 
and only K'^ is zero for LFO, with K<i having a small 
negative value. 







#i 


K[ 


K 2 


KS 


7-8+7- 


LFO 


-158 


7 


-6 







BFO 


-402 


32 









TABLE VIII. Calculated SIA constants (aeV) fitted to Eq.3 
for the Pnma phase of LaFeOs and BiFeOs. 



In Fig. 6 we show a schematic view of the direction of 
the hard axis (red arrows) and easy planes of the 7" 8+ 7" 
structure of LFO. Each local hard axis points in the di- 
rection corresponding to the diagonal of the afr-planes 
in orthorhombic coordinates (#=54.7°, 0=45°). Interest- 
ingly, the hard axis of all Fe sites in successive ac-planes 
point alternately in opposite directions along the b direc- 
tion, as clearly shown in Fig. 6 by the red arrows. We then 
have two possible easy planes (noted A and B on Fig. 6) 
depending of the position of the Fe along the b axis). If 
we look at the combination of the easy planes A and B, 
we find that their intersection is aligned along the c axis. 



It is clear that if we introduce the strong Heisenberg ex- 
change (J ~ 6 meV) that has the tendency to align the 
spins antiparallel, then the intersection of the easy planes 
will determine the spin direction as reported previously 
for Pnma perovskites 30 ' 31 . However, while it has been 
shown from symmetry arguments by Bertaut 12 that the 
SIA allows a canting of the spins in this system, here we 
prove that the combination of all the SIAs compensate 
in such a way that they do not give rise to a canting of 
the spins. 




FIG. 6. Local SIA hard axis (red arrows) of two successive 
ac-planes of Pnma structure (LFO). The intersection of two 
succesive easy planes along the b direction (green and pink 
planes) is indicated by the blue arrow, showing that the two 
different easy planes give rise to an easy axis along the c 
direction. The hard axes point along the [1,1,0] and [1,-1,0] 
directions for the spins in plane B and A respectively. 



E. SIA in a a a R3c and R3c structures 

In Tab. IX we report the SIA parameters obtained for 
BFO and LFO in the R3c (a~a~a~) structure and in the 
R3c structure that includes FE distortions in addition to 
the a~ a~ a~ AFD distortions. As for the Pnma phase, 
in the R3c phase both LFO and BFO have large and 
negative K\ constant (-168 fieV for LFO and -400 /ieV 
for BFO) showing that the SIA is of easy plane form. 
Since the other parameters are equal to zero (K[ and 
K2), we can consider that the SIA in R3c is purely easy 
plane anisotropy. The orientation of this easy plane is 
perpendicular to the hard axis [111]. The spins can thus 
lie freely with any orientation in this easy plane which 
corresponds to the plane of rotation of the experimentally 
observed spin spiral structure of bulk BFO 28 . 

Adding a ferroelectric distortion along the diagonal 
(R3c phase) changes completely the SIA of BFO. The 
amplitude of the anisotropy has the tendency to evolve 
from easy plane [K\ < 0) to easy axis (K\ > 0) as the 
ferroelectric distortions increase. It is interesting to see 
that at the point where the energy is minimum (1.0FE 
in Tab. IX) the SIA is close to the transition where K\ 
changes sign and thus the anisotropy energy is strongly 
reduced. To separate the effect of the FE distortion on 
the SIA, we also performed calculations in the presence of 
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K[ 


K 2 


KS 


9-9-9- 


LFO 


-168 





-1.6 





9-9-9- 


BFO 


-400 





-3.4 





9-9-9-+0.5FE 


BFO 


-281 





- 





9"9"9-+1.0FE 


BFO 


-1.3 










9-9-9-+1.5FE 


BFO 


139 










9-9-9-+1.0FE 


LFO 


-58 










a°a°a° + 1.0FE 


BFO 


217 





-1.9 






TABLE IX. Calculated SIA constants (fieV) obtained by fit- 
ting to Eq.3 of the R3c (9-9-9") and R3c (9"9"9-+FE) 
phases of LFO and BFO. The factor in front of FE describes 
the amplitude of the FE distortions along the [111] direction 
such that the 1.0 amplitude is the point where the energy is 
minimum. For comparison, the last line also shows the SIA 
constants with only the FE distortions (R3m phase), i.e. by 
removing the AFD distortions. 



the FE distortion alone (R3m phase in Tab. IX). The SIA 
is then of easy axis form and it has a large value, close to 
the value of the R3c case (9-9-9-). Therefore, the SIAs 
induced by FE and a~a~a~ distortions are in competi- 
tion in the R3c phase, since a~ a~ a~ favors easy plane 
with K\ < and FE favors easy axis with K\ > 0. The 
resulting SIA in R3c is then determined by the relative 
amplitudes of the FE and AFD distortions. Replacing Bi 
by La in the same structure (9 _ 9~9 - +1.0FE), we also 
find a reduction of the anisotropy but with a less strong 
effect. 

The origin of the differences between La and Bi in the 
SIA in the compounds is difficult to understand. In prin- 
ciple, it could come from the difference of atomic number, 
of cation size or because of the presence of a lone pair on 
the Bi atom. To attempt to isolate these effects, we per- 



formed "computer experiments" by looking at the effect 
of the A-sites on the SIA. We report these in the next 
section. 



F. Differences in SIA between LFO and BFO 

For all the structures for which we calculated the SIA, 
we found the anisotropy to be larger in BFO than in LFO 
by up to a factor of 40. This cannot be attributed to a 
volume effect or to different distortion amplitudes since 
in all cases we adopted the same cell parameters and the 
same amplitude of oxygen rotations or FE distortions. 
This change in the anisotropy can then only originate 
from the differences between Bi and La cations. A first 
difference between Bi 3+ and La 3+ that we can point to is 
the presence of the 6s 2 lone pair on Bi. In principle, this 
could affect the SIA of the magnetic Fe cation since an 
A-site lone pair has the tendency to modify the character 
of the bondlengnth of the B-site cations and anions 34 ' 35 . 
Another important difference is the atomic number (Z) 
which is larger for Bi (83) than for La (57). The ionic 
radii, however, are roughly the same for Bi 3+ and La 3+ 
(around 117 pm). 

In order to understand the A-site effect on the SIA, we 
performed computer experiments in which we replaced 
the A-site in AFe03 by different atoms having a lone pair 
or not and having different atomic numbers and ionic 
radii. For a lone pair candidate we performed calcula- 
tions with Sb at the A-site which has an atomic number 
(51) close to that of La but a smaller ionic radius (90 
pm for Sb 3+ ). For non lone pair candidates, we chose Tl 
which has an atomic number (81) close to that of Bi and 
a small radius (89 pm), and Y which has a much smaller 
atomic number (39) but an intermediate radius (104 pm). 
We report the value of the K\ constant for these different 
A sites and for different structures in Tab.X. 



A cation 


Z 


radius 


lone pair 


a°a°c + 


Ki 

a~ a~ a~ 


a a a +FE 


Y 


39 


104 


no 


16 


-73 


-23 


Sb 


51 


90 


yes 


110 


-161 


12 


La 


57 


117 


no 


7 


-168 


-58 


Tl 


81 


89 


no 


114 


-335 


-57 


Bi 


83 


117 


yes 


264 


-400 


-1 



TABLE X. A atomic number (Z), radius of the A 3+ cation (pm) and SIA constant Ki (/ieV) from Eq.3 of AFe03 with A = 
Bi, La, Tl, Sb and Y within the same frozen structures (a°a°c + , c=10°; a~a~a~ , a=9°; a-a _ a- + 1.0FE, a=9°). The column 
"lone pair" states wether the A 3+ ion has a s lone pair. 



Unfortunately, no obvious trend emerge. In the a°a°c + 
structure, all A-sites have z-easy axis SIA. Tl and Sb 
show similar anisotropy energies, 114 fieV and 110 fieV 



respectively while Y shows a very small anisotropy energy 
of 16 fieV. However, for a~a~a~ there is no obvious trend 
in SIA with atomic number, ionic radius or presence or 
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absence of lone pair, and the origin of the effect of the 
A-site on the SI A is not obvious. We remark that for 
a~ a~ a~ they are yet all easy plane 

If the lone pair has an effect on the SIA, we expect it 
to be more pronounced in the case where the FE distor- 
tions are present since the lone pair is stereochemically 
active in the FE phase. 32 35 . By comparing the last two 
columns in Tab.X, we find indeed that the FE distortion 
causes the largest changes in SIA for compounds with a 
lone pair on the A-site (SbFeOs and BiFeOs). As ex- 
pected from the conclusions in the previous sections, the 
FE distortions clearly decrease the easy plane SIA en- 
ergy of the R3c structure for all the compounds. The 
large SIA energy for BiFeOs m the R3c structure (-400 
fjbeV) is decreased drastically by the FE distortion (-1 
fieV). For SbFeOs the anisotropy is even changed from 
easy plane (-161 /ieV) to easy axis (12 /ieV). Note that 
for consistency in this work we use a U value of 5 eV on 
Fe for all A-site cations. Changes in U value with A-site 
and subsequent effects on the SIA will be the subject of 
future works. 



VI. DISCUSSION 

In the previous sections we have determined the am- 
plitudes of the different magnetic interactions, i.e. ex- 
change, SIA and DM interactions. In the R3c and R3c 
phases, a canting of the spins is allowed by symmetry, 
giving rise to wFM. It has also been demonstrated that 
the only possible mechanism that leads to wFM in R3c 
and R3c is the antisymmetric DM coupling 12 . Our find- 
ing that the calculated SIA gives rise to an easy plane 
perpendicular to the [111] direction and the DM vector 
has the tendency to cant the spins in this easy plane is 
consistent with these earlier symmetry analyses . We 
can conclude that in R3c and R3c, the SIA is neither 
cooperative with nor competing against the DM vector 
since there is almost no SIA energy cost to cant a spin 
in the plane perpendicular to the [111] direction by the 
DM interaction. The direction of canting being the same 
for all the magnetic cations (global wFM), the DM force 
then competes only with the exchange interaction which 
prefers an antiferromagnetic alignment. Thus, the ampli- 
tude of the wFM is directly related to the ratio between 
DM and J such that larger DM and smaller J favor larger 
wFM. 

In the Prima structures, symmetry analysis shows that 
both DM and SIA allow for wFM 12 . Our calculations 
showed that, even if it is allowed by symmetry, the SIA 
does not give rise to wFM. While the local SIA of different 
spin sublattices are not the same, the sum of the SIAs on 
all the spin sublattices gives rise to a fixed unidirectional 
orientation of the spins (global easy axis or MCA) . As we 
saw in Fig. 6, the local easy planes are not parallel to the 
direction where the DM forces want to cant the spins. 
This means that in the Prima structure, the canting of 
the spins is, as in the R3c structure, only due to the 



DM interaction. In contrast to the R3c case, where the 
SIA does not affect the development of the wFM, in the 
Prima structure the SIA competes with the DM. As a 
result, in the Prima phase, the DM has to compete with 
not only the exchange interaction but also the SIA. 

Our results allow us to conclude that depending on 
the structure, the amplitude of the canting of the spins 
in Prima and R3c phases comes from a delicate balance 
between SIA, DM and J. If one wants to design large 
wFM through magnetic interactions in these distorted 
perovskites, one has to reduce J and SIA relative to the 
DM interaction. The main parameters which we can play 
with to design large canting are the amplitude of the FE 
and AFD distortions of the perovskite structure. To that 
end, we need to understand the coupling between the 
crystal distortions and the canting of the spins, a link we 
analyze in the next section. 



VII. EFFECT OF AFD AMPLITUDE ON J, DM 
AND SIA 

As seen in the previous sections, it is clear that AFDs 
have the tendency to reduce the exchange parameter J, 
to induce large SIA and to allow for DM in the perovskite 
structure. As a result, they give rise to wFM in Prima 
and R3c phases for example. In our calculations above, 
however, we calculated the amplitude of the wFM for 
only one amplitude of AFD distortions in each crystal 
phase. Here, we investigate how the amplitude of the 
crystal distortions acts on the non-collinear magnetism 
and thus on the wFM. We thus analyze the influence 
of the distortions on the non-collinear magnetism in the 
Prima and R3c phases by calculating the evolution of 
J, DM, SIA and wFM as a function of the amplitude of 
AFD and FE distortions. 

In Fig. 7 (a) we report the evolution of the spin cant- 
ing angle responsible for the wFM as a function of the 
AFD amplitudes in LFO in the R3c phase. For 0° to 
about 10° of AFD distortion, the wFM moment increases 
linearly with increasing rotation (red curve). In Fig. 7 
(b), (c) and (d) we report the evolution of respectively 
the total exchange parameter J, the SIA energy between 
the hard and easy direction and the total DM vector of 
the R3c structure as a function of the oxygen octahedral 
angle of rotation. The total J is just the sum over all 
six nearest-neighbor exchange constants and D Tot is the 
norm of the total DM vector. As expected, J decreases 
with the increase of the amplitude of the oxygen octahe- 
dra rotation, and we find that the decrease is linear for 
rotation angles from 0° to 15°. In the same range, the 
DM increases linearly with the rotation angle and the 
SIA increases quadratically. In Fig. 7 (a) we also report 
the value of arctan(Z)/ J) which Interestingly follows per- 
fectly the amplitude of the wFM. This means that, when 
the SIA does not compete with DM, the amplitude of the 
spin canting can be directly related to the ratio between 
DM and J and thus to the oxygen rotation amplitude. 
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FIG. 7. Evolution of (a) the canting angle responsible for 
the wFM (red circles) (b) the exchange parameter J, (c) the 
SIA energy difference between hard and esasy direction of 
the spin and (d) the amplitude of the DM vector versus the 
oxygen rotation angle in i?3c-structure LFO. In plot (a) we 
also report the values of aictaxi(DTot/ Jrot) (blue triangles). 



In Fig. 8 (a) we report the evolution of the spin canting 
angle responsible for the wFM versus the AFD ampli- 
tudes of LFO in the Prima phase. Fixing the out-of- 
phase rotations and varying only the amplitude of the 
in-phase rotation causes almost no change in the am- 
plitude of the wFM (blue squares in Fig. 8 (a)). How- 
ever, freezing the in-phase rotations and varying the out- 
of-phase rotations changes the magnitude of the wFM, 
clearly showing that the wFM is directly linked to the 
amplitude of the out-of-phase rotations (green circles in 
Fig. 8 a). For the whole range of out-of-phase rotation 
angles considered (from 0° to 15°), we observe that the 
relationship between AFD and wFM is linear while it de- 
viates from perfect linear behaviour in the R3c structure 
beyond oxygen rotation angle of about 12° (Fig. 7 (a)). 

In Fig. 8 (b) and (c) we report the evolution of SIA and 
the component of the DM vector that is responsible for 
wFM (the x-component as reported in section. IV A). As 
can be seen from the blue squares in Fig. 8 (b) and (c), 
the in-phase AFD has almost no effect on SIA and the 
DM interaction. Only the out-of-phase oxygen rotations 
have a sizeable effect on the SIA and the DM interaction 
where, as for i?3c, SIA and D^ot increase with increas- 
ing out-of-phase AFD distortions. This suggests that, to 
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FIG. 8. Evolution of (a) the canting angle responsible for the 
wFM (b) the SIA energy difference between hard and esasy 
direction of the spin and (c) the amplitude of the DM vector 
versus the oxygen rotation angle of the in-phase rotation (blue 
squares) and the two out-of-phase rotations (green circles) of 
the a~b^~a~ structure of LFO. When the out-of-phase angle is 
changes, the in-phase angle is fixed at 8° (x~8 + x~) and when 
the in-phase rotation is varied, the out-of-phase rotations are 
fixed at 7° (7"x + 7"). 



VIII. CONCLUSION 

In this paper we analyzed the SIA, DM and exchange 
interactions in perovskite structures BiFeOa an d LaFeOa 
with a range of structural distortions. We analyzed the 
effect of these three interactions on the final magnetic 
ground state in different distortion patterns that can 
be present in the perovskites: a°a°c + , a°a°c _ , a°6 + 6 + , 
a°b~b~ , a~6 + a~, a~ a~ a~ and FE distortions. We con- 
firmed that in all cases the spin canting is due only to 
the DM interaction. This is true even in the a~b + a~ 
case where symmetry allows also a possible contribution 
from the SIA 12 . We found that the amplitude of the 
canting of the spins, and thus of the wFM, is determined 
by a balance between the amplitude of the DM inter- 
action, the SIA and the exchange interaction. The ex- 
change interaction is always in competition with DM for 
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the Fe 3+ of the studied compounds, while the SIA can 
be neutral, cooperative or in competition with DM. For 
both chemistries we found that the SIA is not impor- 
tant in the non-ferroelectric R3c and Prima structures. 
This is in contrast to the case of Prima LaMnOa, where 
SIA seems to be the important mechanism that leads to 
wFM 36 . 

We also looked at how the amplitude of the atomic 
distortions affect the different magnetic interactions and 
found that the change of the DM and the J's are close 
to linear in the AFD amplitudes, while the SIA changes 
more drastically with these amplitudes. The J values 
are reduced by the AFD distortions while the DM vec- 
tors and SIA increase with the AFD distortions. The 
amplitude of the wFM is then directly linked to the ratio 
between the amplitude of the DM interaction and the J 
values such that the wFM canting angle follows a geo- 
metrical law with respect to this ratio and the angle of 
the oxygen rotations in the R3c and Prima cases. It was 
also interesting to see that in the Prima case, mixing in- 
phase and out-of-phase AFD (a~6 + a~), the wFM is only 
induced by the out-of-phase AFD. Thus, to control the 
wFM in the Prima type of structures, one has to tune 
the out-of-phase AFD. 

By decomposing the magnetic interactions in BFO and 
LFO, we can conclude that the DM and J parameters can 
be simply linked to the atomic distortions while the SIA 
has a much more complex behavior with respect to lattice 
modifications and is much more sensitive to them. We 
note that the SIA values we computed here were all for 
the d 5 configuration of Fe 3+ . While a simple local atomic 
environment view of the SIA and crystal field splitting 
would lead us to predict small anisotropics (L = and 
thus small spin-orbit interaction) for d 5 electrons in "cu- 
bic" local oxygen octahedra (the Fe-0 distances were 
kept constant and equal in our calculations) , surprisingly 
we found for some structures very large anisotropics of 
the order of several hundred /ieV. This was mainly re- 
lated to an A-site effect in the ABO3 structure. The 
size, the atomic number and the presence of a lone pair 
on the A-site seem to have a strong influence on the SIA 
and are responsible for the large anisotropics reported 
for the Fe 3+ . Since in our calculations we found a very 
small value of the orbital magnetic moment (<0.02 lib), 



and since the the spin-orbit interaction term has the form 
XL • 5, the large induced SIA from the A cation must be 
through the change in the spin-orbit coupling constant A. 
More complex interplay between the lattice distortions, 
the chemistry of the cations and the SIA can be expected 
when going away from d 5 electronic configuration of the 
B cation since the SIA will also originate from aL/0 
contribution in the spin-orbit interaction 37,38 . 

Finally, we found that FE distortions, as present in the 
ground state of BFO, have a smaller effect on the DM and 
J constants and thus on the wFM. However, we found 
that FE distortion has a strong influence on SIA since FE 
distortion favors an easy axis type of SIA which is in com- 
petition with the easy plane anisotropy induced by AFD 
distortions in the R3c structure. This could be particu- 
larly interesting for technological applications for the fol- 
lowing reason: the FE distortions tend to orient the spins 
along the [111] direction and the AFD distortions in the 
plane perpendicular to the [111] direction. Then depend- 
ing of the balance between the amplitude of the AFD and 
the amplitude of the FE (AFD/FE), the spins can lie ei- 
ther in the plane perpendicular to the [111] direction or 
parallel to the [111] direction. When the spins lie in the 
[111] direction, no wFM is allowed by symmetry (because 
the DM vector is also along the [111] direction 13 ). Then, 
if one can simply tune the ratio AFD/FE in BFO by 
acting on the amplitude of these distortions (for exam- 
ple with pressure 39 , epitaxial strain 40 42 , electric field 43 
or phononic exitation), the wFM could be switched on 
and off through the competition between the SIA shapes 
(easy-plane versus easy-axis) of the AFD and FE distor- 
tions. This process could be more readily attainable than 
the reversal of wFM that has been discussed previously 6 ' 8 
and opens the door to new design rules for controlling 
magnetism by nonmagnetic external parameters. 
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